EULERIAN PARTITIONS FOR CONFIGURATIONS OF 

SKEW LINES 



ROLAND BACHERi AND DAVID CAREER^ 

Abstract. In this paper, which is a complement of [I], we study a 
few elementary invariants for configurations of skew lines, as intro- 
duced and analyzed first by Viro and his collaborators. We slightly 
simplify the exposition of some known invariants and use them to 
define a natural partition of the lines in a skew configuration. 

We also describe an algorithm which constructs a spindle-permutation 
for a given switching class, or proves non-existence of such a spindle- 
permutation. 



1. Introduction 

A configuration of n skew lines in or a skew configuration or 
an interlacing of skew lines is a set of n non-intersecting lines in 
containing no pair of parallel lines. 

Skew configurations are only considered up to rigid isotopies (con- 
tinuous deformations of skew configurations or, equivalently, isotopies 
of the ambient space under which lines remain pairwise skew lines). 

The study and classification of configurations of skew lines up to 
isotopy was initiated by Viro [12] and pursued by Viro, Mazurovskii, 
Borobia-Mazurovskii, Drobotukhina and Khashin, see for example [21 
131 [5l El [HI m [13] . The survey paper [13] (and its updates available on 
the authors web-sites and on the arXiv) contains historical information 
and is a good introduction into the subject and its higher-dimensional 
generalizations. Most of these results are also exposed in the survey 
paper [3] from which we borrow some terminology not used in the 
original work of Viro's school. 

A spindle (or isotopy join or horizontal configuration) is a particu- 
larly nice configuration of skew lines in which all lines intersect an ori- 
ented additional line A, called the axis of the spindle. Its isotopy class 
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is completely described by a spindle-permutation a : {l,...,n} — 
{1,...,?T,} encoding the order in which an open half-plane revolving 
around its boundary A intersects the lines during a half-turn (see Sec- 
tion [5] for the precise definition). A spindle- configuration is a skew 
configuration isotopic to a spindle. 

Three types of combinatorial moves (described in Section [5]) of a 
spindle permutation yield isotopic spindles-configurations and gener- 
ate an equivalence relation, called the spindle- equivalence relation, on 
permutations of {1, . . . , n}. 

Isotopy classes of spindle-configurations are well understood by the 
following combinatorial description, given in pQ: 

Theorem 1.1. Two spindle-permutations a, a' give rise to isotopic 
spindle-configurations if and only if a and a' are spindle-equivalent. 

Orienting and labeling all lines of a skew configuration, one gets 
a linking matrix encoding isotopy classes for pairs of oriented skew 
lines. The associated switching class or homology equivalence class 
is independent of labels and orientations. A result of Khashin and 
Mazurovskii [Tij Theorem 3.2] states that homology-equivalent spindles 
(spindles defining the same switching class) are isotopic. Thus we have: 

Corollary 1.2. Two spindle-permutations define the same switching 
class if and only if they are spindle- equivalent. 

In this paper we define the Euler partition for a switching class. The 
definition depends on the parity of the order and can be refined to an 
Euler tree for switching classes of even order. 

The last part of the paper describes an algorithm for computing a 
spindle-permutation (or proving its non-existence) for a given switching 
class. 

The sequel of this paper is organized as follows: Sections [2] and [3] are 
devoted to (various aspects of) switching classes. Section H] deals with 
the Eulerian partition induced by the switching classes. Section O intro- 
duces spindle-structures. In Section E] we describe an algorithm which 
computes (or proves non-existence of) a spindle-permutation (which 
is unique up to spindle-equivalence by Corollary 11.21) having a linking 
matrix of given switching class. 

2. LiNKINC MATRICES AND SWITCHINC CLASSES 

Pairs of oriented under- or over-crossing curves (as arising for in- 
stance from oriented knots and links) can be encoded by signs ±1 as 
shown in Figure [H 
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positive 
crossing 



negative 
crossing 



Figure 1. Positive and negative crossings 



The sign or linking number lk{L a, Lb) G {=tl} between two oriented 
skew lines La,Lb C was introduced by Viro [12]. The linking 
matrix of a configuration involving n oriented and labeled skew lines 
Li, . . . ,Ln is the symmetric n x n matrix X with diagonal coefficients 
Xi^i = and Xij = lk(Lj, Lj) G {±1} for i ^ j. 



Figure 2. A configuration of 6 labeled and oriented 
skew lines 

Figure [2] shows a labeled and oriented configuration of six skew lines 
with linking matrix 
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Two symmetric matrices X and Y are switching-equivalent if 

Y = D X P D 

where P is a permutation matrix and D is a diagonal matrix with rfj j G 
{±1}. Since P and D are orthogonal, we have {PD)^^ = D^P* = DP*. 
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Switching-equivalent matrices are thus conjugate and have the same 
characteristic polynomial. 

Proposition 2.1. All linking matrices of a fixed configuration of skew 
lines are switching-equivalent. 

Proof. Relabeling the lines conjugates a linking matrix X by a per- 
mutation matrix. Inverting the orientation of some lines amounts to 
conjugation by a diagonal ±1 matrix. □ 

Remark 2.2. The terminology "switching classes" (many authors speak 
of ''two-graphs" which is the standard terminology for the underlying 
combinatorial object) is motivated by the following combinatorial inter- 
pretation and definition of switching classes. 

Two finite simple (loopless and no multiple edges) graphs Ti = (V, Ei) 
and r2 = (V, E2) are switching-related with respect to a subset of ver- 
tices V^^ C V if their edge-sets Ei,E2 coincide on {y^xVJ)U{{y\VJ)x 
{y\V-)') and are complementary on (V^_ x {y\V-)) U {{y\V-) x V-) . 
A switching class of graphs is an equivalence class of switching-related 
graphs, see Figure\^for two graphs in a common switching class. 

12 12 



34 34 

FicURE 3. Fi and F2 are switching-related with respect 
to {1,2} C {1,2,3,4} 

Encoding adjacency, respectively non- adjacency, of distinct vertices 
by 1, respectively —1, yields a bijection between switching classes of 
graphs and switching classes of matrices. Conjugation by permutation 
matrices corresponds to relabeling the vertices of a graph F. Conjuga- 
tion by a diagonal -^1 — matrix corresponds to the substitution of T by 
a switching-related graph. 

3. Switching classes and vorticity 

The set of vorticities (also called homological equivalence class or 
chiral signature), introduced by Viro [12], is a classical and well-known 
invariant for configurations of skew lines. We sketch below briefly the 
well-known proof that it corresponds to the switching class of an asso- 
ciated linking matrix. 
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We prefer to work with switching classes corresponding to symmetric 
matrices (up to conjugation by signed permutation matrices) with zero 
diagonal and off-diagonal coefficients in {±1}. 

The vorticity vort(Lj, Lj, L^) of three lines (see ^2] Section 2] or [H 
Section 3]) is defined as the product XijXj^kXk,i G {±1} of the signs 
for the corresponding three crossings. The result is independent of 
the chosen orientations for Li,Lj,Lk, classifies the skew-configuration 
{Li, Lj, Lk} up to isotopy and yields an invariant 

{triplets of lines in configurations of skew lines} — > {=tl} • 

Let us remark that almost all authors use the terminology linking 
coefficient instead of vorticity. This is slightly confusing since the link- 
ing coefficient denotes also the isotopy type of a pair of oriented skew 
lines. 

The set of vorticities is the list of vorticities vort(Lj, L^-, L^) for all 
triplets of lines {Li, Lj, L^] in a configuration of skew lines. 

Sets of vorticities (defining a two- graph, see [H]) and switching classes 
are equivalent. Indeed, vorticities of a configuration of skew lines C can 
easily be retrieved from a linking matrix for C. Conversely, given all 
vorticites vort(Lj, Lj, L^) of a configuration C = {Li, . . . , Ln}, choose 
an orientation of the first line Li. Orient the remaining lines L2, . . . ,Ln 
such that they cross Li positively. A linking matrix X for C is given 
by Xi^i = Xi^i = 1 for i such that 2 < i < n and Xa^b = vort(Li, La, Lh) 
for 2 < a ^ b < n. 

Two configurations of skew lines are homologically equivalent (the 
terminology refers to properties of the complement, endowed with a 
suitable extra-structure, of a configuration in MP^) if there exists a bi- 
jection between their lines, which preserves all vorticities. Two config- 
urations are homologically equivalent if and only if they have switching- 
equivalent linking matrices. 

The sign indeterminacy in linking matrices representing switching 
classes makes their use more difficult. A satisfactory answer addressing 
this problem will be given in Section 14.11 for switching classes of odd 
order. For even orders, there seems to be no completely satisfactory 
way to get rid of all sign-indeterminacies, see Section 14.21 

4. EULER PARTITIONS 

In this section we study invariants of computational cost O(n^) for 
switching classes of order n. 

The behaviour of switching classes depends on the parity of their 
order. 
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Switching classes of odd order 2n — 1 are in bijection with Eulerian 
graphs. This endows the hnes of a skew configuration consisting of an 
odd number of hnes with a semi- orientation (a canonical orientation, 
up to global change) which we call the Eulerian semi-orientation. An 
Eulerian semi-orientation induces a partition of the lines into equiva- 
lence classes by counting their number of positive crossings. We con- 
sider the case of odd order in Section 14. 1[ 

The situation for switching classes of even order 2n is more com- 
plicated. We replace Eulerian graphs appearing for odd orders by a 
suitable kind of rooted binary trees which we call Euler trees. The 
leaves of the Euler tree induce again a natural partition, called the 
Euler partition, of the set of lines into equivalence classes of even car- 
dinalities. Section 14.21 deals with the even case. 

4.1. Switching classes of odd order - Eulerian semi-orientations. 

A simple finite graph F is Eulerian if all its vertices are of even degree. 
Figure H] shows all seven Eulerian graphs on 5 vertices. 



The following well-known result goes back to Seidel [lU] . 

Proposition 4.1. Eulerian graphs with an odd number 2n — 1 of ver- 
tices are in bijection with switching classes of order 2n — 1. 

We recall a simple proof of Proposition 14.11 since it yields a fast 
algorithm for computing Eulerian semi-orientations on configurations 
with an odd number of skew lines. 

Proof. Choose a representing matrix X of a switching class. For i such 
that 1 < i < 2n — 1 define the number 




Figure 4. All Eulerian graphs on 5 vertices 




2n-l 
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counting all entries equal to 1 in the i—th row of X. Since X is sym- 
metric, the vector {vi, . . . , V2n-i) has an even number of odd coefficients 
and conjugation of the matrix X with the diagonal matrix having di- 
agonal entries (—1)^' turns X into a matrix Xe having an even number 
of I's in each row and column. The matrix Xe is well-defined up to 
conjugation by a permutation matrix and defines an Eulerian graph 
Fx with vertices {1, . . . , 2n — 1} and edges {i,j} if {XE)i,j = 1- The 
Eulerian graph Fx is unique up to relabeling its vertices since switch- 
ing with respect to a non-trivial subset of vertices destroys the Eulerian 
property of Tx- □ 

A semi- orientation of a set of lines C is an orientation of all lines in 
£, up to global inversion of all orientiations. 

Let C be a configuration of skew lines having an odd number of lines. 
Label and orient the lines of C arbitrarily in order to get a linking 
matrix X. Inverting the orientations of all lines having an odd number 
of positive crossings we get a unique semi-orientation which we call the 
Eulerian semi- orientation of C. 

An Eulerian linking matrix Xe associated to an Eulerian semi-orientation 
of C is uniquely defined up to conjugation by a permutation matrix. 
Its invariants coincide with those of the switching class of Xe but are 
slightly easier to compute since there is no sign ambiguity. In particu- 
lar, some of them can be computed using only O(n^) operations. 

An Eulerian partition of the set of lines of a configuration consisting 
of an odd number of lines is by definition the partition of the lines 
into subsets Ck consisting of all lines involved in exactly 2k positive 
crossings for an Eulerian semi-orientation. 

A few more invariants of Eulerian matrices are: 

(1) The total sum j of all entries in an Eulerian linking ma- 
trix Xe (this is of course equivalent to the computation of the 
number of entries equal to 1 in Xe)- The computation of this 
invariant needs only O(n^) operations. 

(2) Its signature e = Y[i<j^i,j- The easy identity 



relates the signature to the total sum ■ Xij of all entries in 
an Eulerian linking matrix. 
(3) The number of rows of Xe with given row-sum. These numbers 
yield of course the cardinalities of the sets Co,Ci, . . . and can 
be computed using O(n^) operations. 
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(4) All invariants of tlic associated Eulerian graph (having edges 
corresponding to entries Xij — 1) defined by Xe, e.g. the num- 
ber of triangles or of other fixed subgraphs. In particular, one 
can consider the number ^ of edges joining a vertex of degree 
2i to a vertex of degree 2j. 

For example, for 7 vertices, there are 54 different Eulerian graphs, 
36 different sequences of vertex degrees (up to a permutation of the 
vertices), and 18 different numbers for the cardinality of I's in Xe- 

4.2. Svi^itching classes of even order - Euler partitions. The sit- 
uation in this case is more complicated and less satisfactory. 

Given a matrix X representing a switching class with 2n vertices, 
there exists a natural partition of the 2n rows it! of X into two subsets 
R+ and R- according to the sign 

associated to the i—th row of X. This sign is well-defined since switch- 
ing (conjugation) with respect to a diagonal {±1}— matrix D multiplies 
both factors Ylj^i^ij Y[s<t-'^s,t by det(i5) G {±1}- 

Since Hi = Ui,j, i^j ^ij {Ils<t^s,tY'' = 1, both subsets R+,R- 
have even cardinalities. 

If (or equivalently, i?_) is non-empty, it defines a symmetric 
submatrix of even size ^{R+) corresponding to all rows and cohimns 
with indices in Iterating the above construction we get a partition 
i?+ = U This construction is most conveniently encoded by 

a rooted binary tree embedded in the oriented plane which we call the 
Euler tree of X: Draw a root R corresponding to the row-set R oi X. 
If the partition R = R^ U R is non-trivial, join the root i? to a left 
successor called i?_ and a right successor called The Euler tree of 
X is now constructed recursively by gluing the root R± of the Euler 
tree associated to X± onto the corresponding successor R± of the root 
R. 

The leaves of the Euler tree T{X) of X correspond to subsets R^ 
(with w G {±}*) of even cardinalities 2n^ summing up to 2n. Leaves 
of T{X) define symmetric submatrices in X which we call Eulerian: 
All their row-sums are identical modulo 2 and can be chosen to be 
even, after a suitable conjugation. The row partition R = i?+ U i?_ 
of an Eulerian matrix X of even size is by definition trivial. The sign 
e G {±1} defined as e = for i an arbitrary row of X is called the 
signature of the Eulerian matrix X. An Eulerian matrix of size 2 
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has always signature 1. For Eulerian matrices of size 2n > 4 both 
signs can occur as signatures since changing the signs of the entries 
Xij, 1 < "j 7^ J < 3, inverts the signature of an Eulerian matrix (and 
preserves the set of Eulerian matrices of even size > 4). The signature 
of an Eulerian matrix encodes the parity of the number of edges in an 
Eulerian graph (having only vertices of even degrees) in the switching 
class of X. 

The leaves of the Euler tree define a natural partition of the set of 
rows of X into subsets. We call this partition the Euler partition. 

An Euler tree is signed if its leaves are endowed with signs ±1 cor- 
responding to the signs of the associated Eulerian matrices. An Euler 
tree is weighted if its leaves are endowed with strictly positive natu- 
ral weights, a weight m corresponding to an Eulerian matrix of size 
2m X 2m. A signed weighted Euler tree is both signed and weighted. 

Example 4.2. The symmetric matrix 
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yields the Euler partition 

i?__ = {3,5} U i?_+ = {6,10} U i?+ = {1,2,4,7,8,9} 

(where the Eulerian submatrix associated to i?+ has signature 1). The 
associated signed Euler tree (with leaves of respective weights 1, 1 and 
3) is depicted in Figure O 




Figure 5. The signed Euler tree defined by Example 14.21 
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Example 4.3. The characteristic polynomial of a linking matrix of a 
configuration of skew lines is in general weaker than its switching class: 
The linking matrices 
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are in different switching classes: The row-partition R = i?_ U i?+ 
of the first matrix is given by R- = {2,5,7,8} and i?+ = {1,3,4,6} 
with associated Eulerian matrices and X+ both of signature 1 . The 
second matrix is Eulerian with signature 1. On the other hand, they 
have the same characteristic polynomial 



{t -3){t- lf{t + l){t + 3)^(i^ -2t- 11) . 



This example is minimal in the sense that distinct switching classes of 
order less than 8 have distinct characteristic polynomials. 

Let us mention a last invariant related to the Euler tree for a switch- 
ing class X having even order 2n. Let Ri, . . . , Rm <Z Rhe the Euler 
partition of X . For 1 < i < j < m, define numbers a^j e {±1} by 



where Si^ G i?i is a fixed element. Note that the number ai^i is the 
signature of the Eulerian matrix defined by the rows (and columns) of 
the set Ri- One can easily check that the numbers well-defined. 



and 
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Remark 4.4. The equivalence relation induced on lines by the Euler 
partition is fairly coarse. It is for instance generally much rougher that 
the equivalence relation given by homologous lines defined by Viro [12j. 

4.3. Enumerative aspects. It is natural to enumerate (signed) weighted 
Euler trees according to the total sum n of all weights. 

The generating function F{z) = Yl'^=o ^nZ'^ enumerating the number 
an of distinct weighted Euler trees with total weight n satisfies the 
equation 

F{z) = + {F{z) - If 

(with ao = 1 corresponding to the empty tree). Indeed, weighted 
Euler trees reduced to a leaf contribute 1/(1 — z) to F{z). All other 
weighted Euler trees are obtained by gluing two weighted Euler trees of 
strictly positive weights below a root and are enumerated by the factor 
{F{z) - If. 

Solving for F{z) we get the closed form 

:i{l-z)-^{l-z){l-hz) 



m = E 



n=0 



""^ 2(1 -z) 



showing that 



hm = 5 . 



n->-oo an 

The first terms ao,ai, . . . are given by 

1, 1, 2, 5, 15, 51, 188, 731, 2950, 12235, . . . , 

see also Sequence A7317 of [TTj. 

Similarly, the generating function Fs{z) = ^^o-^"-^" enumerating 
the number /3„ of signed weighted Euler trees (keeping also track of the 
signature of all leaves with weight > 2) with total weight n satisfies 
the equation 

v2 



Fsiz) = ^ + (F^iz) - If 



We get thus 



n=0 



and 



lim %1 = ~ 5.7016 



The first terms /3o, • • • are given by 



1, 1, 3, 8, 27, 104, 436, 1930, 8871, 41916, . . . , 
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see Sequence A110886 of [II]. 

5. Spindle structures for switchinc classes 

A construction of Mazurovskii originally called the isotopy join (or 
spindle) associates a configuration of n skew lines to every permutation 
of n letters. We recall that a spindle is a configuration of skew lines 
with all lines intersecting an oriented auxiliary line A, called its axis. A 
spindle- configuration (or a spindle structure) is a configuration of skew 
lines isotopic to a spindle. 

The orientation of the axis A induces a linear order Li < ■ ■ ■ < L„ 
on the n lines of a spindle C. Each line & C defines a plane Ilj 
containing Lj and the axis A. 

A second oriented auxiliary line B (called a directrix) in general posi- 
tion with respect to A, IIi, . . . , n„ and crossing A negatively, intersects 
the planes Hi, . . . ,Iln at points a{Li) = B n Hi. One can assume 
a{Li) G by a suitable rotation fixing ACiUi of the plane flj contain- 
ing Li. Since the orientation of B induces a linear order on the points 
(T{Li), we get a spindle-permutation (still denoted) i \ — )■ a{i) of the set 
{1, . . . ,n} by identifying the two hnearly ordered sets Li, . . . , L„ and 
cr(Li), . . . , a{Ln) in the obvious way with {!,..., n}. Figure E] displays 
an example corresponding to cr(l) = 1, a (2) =4, a (3) = 2, a (4) = 
5, (t(5) = 3. 



A linking matrix X of a spindle C is easily computed as follows. 
Transform C isotopically into a spindle with oriented axis A and di- 
rectrix B as above. Orient a line Lj from Lid A to (j{Li) = LiCi B. A 
straightforward computation shows that the linking matrix X of this 
labeled and oriented configuration of skew lines has coefficients 



A 




Figure 6. A spindle 



Xij = sign((i - j){(j{i) - a{j))) 
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where sign(O) = and sign(a;) = -j^ for x 7^ and where a is the 
corresponding spindle-permutation. 
The hnking matrix of Figure |6] is 



X 



Two spindle-permutations are equivalent (see jH Section 15]), and 
give rise to isotopic spindle configurations, if they are equivalent under 
the equivalence relation generated by 

(1) (Circular move) 

0" ~ yU if /i(z) = (s + cr((z + t) (mod n))) (mod n) 

for some integers < s, t < n (all integers are modulo n). 

(2) (Vertical reflection of a block or local reversal) a ~ /i if cr([l, k]) = 
[1, k] and 



fi{i) 



k + 1 - a{k + I 



(Tit 



i < k 
i > k 



for some integer k <n (see Figure [7]). 







123456 123456 

Figure 7. Vertical reflection of a block 

(3) (Horizontal reflection of a block or local inversion) a ~ yU if 
there exists an integer 1 < k < n such that cr([l,A;]) = [1,A;] 
and 

a-^{i) i < k 

a(i) i > k 



(see Figure [HD- 

Permutations giving rise to linking matrices in a common switching 
class are spindle-equivalent by Corollary 11.21 
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123456 1234 





123456 12345 

Figure 8. Horizontal reflection of a block 

6. An ALGORITHM FOR A SPINDLE-STRUCTURE 

Given a switching class represented by a matrix X, the following 
algorithm constructs a spindle-permutation (which is unique up to 
spindle-equivalence) with linking matrix in the switching class of X 
or proves non-existence of such a permutation. 

Algorithm 6.1. 

Initial data. A natural number n and a switching class represented by a sym- 
metric matrix X of order n with rows and columns indexed by 
{1, . . . , n} and coefficients Xij satisfying 



X. 



0, l<i<n , 

Xj^i e {±1} , I <iy^ j <n . 



Initialization. Conjugate the symmetric matrix X by the diagonal matrix with 
diagonal coefficients {l,Xi^2,Xi,3, ■ ■ ■ ,Xi^n)- Set = 7(2) = 
1, a(l) = 1 andk = 2. 
Main loop. Replace ^{k) by ^{k) + 1 and set 

k-l 

a{k) = 1 + tl{j I x^(fc),j = -1} + ^x^(s),'y{k) ■ 

s=l 

Check the following conditions: 

(1) j{k)j^^{s)forse{l,...,k~l}. 

(2) x^(k),-,is) = sign{a{k) - a{s)) /or s G {1, . . . , /c - 1} (where 
sign{0) = and sign{x) = |^ for x ^ 0). 

(3) forj e {l,...,n}\{7(l),...,7(A;)} and for s e {1,...,A;- 

^/ •^j,7(s) *^7(s),7(fc) 1; then 3^_;,7(fe) — 

If all conditions are fulfilled, then: 

if k = n, print all the data (mainly the spindle-permutation 
i I — > a{i) and perhaps also the conjugating permutation 
i I — > 7(i) ) and stop. 
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if k < n, then set 'j{k + 1) = 1, replace k by k + 1 and 
iterate the main loop. 

If at least one of the above conditions is not fulfilled, then: 
while 7 (A;) = n replace k by k — 1. 

if k = 1, print "no spindle structure exists for this switch- 
ing class" and stop, 
if k > 1, iterate the main loop. 

6.1. Explanation of the algorithm. The initialization is simply a 
normalization: we assume that the first row of the matrix represents 
the first line of a spindle-permutation a normalized to a{l) = 1 (up to 
a circular move, this can always be done for a spindle-permutation). 

The main loop assumes that row number '^{k) of X contains the 
linking numbers of the A;— th line Lk (supposing a correct possible choice 
of the rows encoding the linking numbers of Li, . . . , The image 

cr(/c) of k under a spindle-permutation is then uniquely defined and 
given by the formula used in the main loop. 

One has to check three necessary conditions: 

• The first condition checks that row number 7(/c) has not been 
used before. 

• The second condition checks the consistency of the choice for 

'j{k) with all previous choices. 

• If the third condition is violated, the choice of rows 7(1), . . . , ■j{k) 
leads to a dead end. Indeed, we have then either = 
l,Xj^j(^s) — ~1 or x-y(5)^-y(fc) = —l,Xj^^(^s) — 1 for some index 
j ^ {7(1), 7(2), . . . ,7(A;)} and some natural integer s < k. In 
the first case, the line-segments Lg and Lk with s < k represent- 
ing a graphically do not cross. This shows that any line-segment 
Lm with m > k which crosses Lg has to cross Lk first. Since 
there must be at least one such line segment corresponding to 
the choise 7(7) = m, the algorithm must backtrack. The second 
case is similar. 

The algorithm runs correctly even without checking out Condition 
(3). However, it loses much of its interest: An instance of Condition 
(3) (with fixed j, s, k) is violated with probability | for a "random" 
choice (made e.g. by hipping a fair coin) of Xj^^(^s)-iXj,'y{k) £ {il}- This 
ensures fast running time in the average, as observed experimentally. 

The algorithm, if successful, produces two permutations a and 7. 
The linking matrix of the spindle permutation cr is in the switching 
class of X and 7 yields a conjugation between these two matrices. 
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More precisely: 

^i(i)Mj) = sign((i - - o^(j))) 

under the assumption xi^i = Xi^i = 1 for 2 < i < n. 

Failure of the algorithm (indicated by the output "no spindle struc- 
ture exists for this switching class") proves non-existence of a spindle 
structure in the switching class of X. 
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